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I. INTRODUCTION 



There is a class of electroweak contributions to the muon g ~ 2 containing a fermion 
triangle along with a virtual photon and Z boson as shown at Fig. ^that is discussed in 
I^M 1^; IE ly| ■ For the determination of the muon anomalous electromagnetic moment we 




FIG. 1: Effective Zjj* coupling contributing to aj^^ ■ 

are interested in the Z* 7* transition in the presence of the external magnetic field to 
first order in this field. In this approximation one may consider the current as a source 
of a soft photon with polarization vector e^{k) and momentum /c — *• 0. In such kinematics 
projection of the amplitude on e'^{k) contains only two Lorentz- invariant structures 

+ WL{p^,rr?)po,pPfp^, (1.1) 

f pLV "^^fiwysf ) f/iu kpCy kyCp^ 

and can be viewed as a correlator of the axial and vector currents in the external electro- 
magnetic field with strength tensor f^j^. The same expression appears while analyzing the 
dominant contribution of light-by-light scattering to g — 2. 

Both structures wt and wl are transversal with respect to vector current, but only the 
first structure is transversal with respect to axial current while the second is longitudinal. 
According to the classical papers by Rosenberg iS], Adler 0], Bell and Jackiw 9] at the 
one-loop level the invariant functions wl^t satisfy the relation 

= 2w^^-'"°''\p\7n'). (1.2) 

It was shown by A.I. Vainshtein that there is the symmetry of the triangle amplitude 
Taixv under permutation z/ ^ a in the chiral limit. This symmetry preserves the relation 
between wl and wt (!1.2p for the case m = in any order of perturbation theory. 

Moreover, according to the Adler-Bardeen theorem jllj the anomalous longitudinal part 
of the triangle is not renormalized in the chiral limit. It is worthy to note that this state- 
ment implies an operator relation, while the matrix elements get the corrections l3| due 



to anomalous dimension of axial current. At the same time, the validity of Adler-Bardeen 
theorem at the operator level allows to express jl3| these three-loop corrections ^ in terms 
of earlier two-loop and even one-loop calculations of anomalous divergencies. 
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To apply Adler-Bardeen theorem to the problem in question one should recall that the 
axial anomaly is expressed only through the longitudinal part |1, 0| 

p"T,,(/,m2 = 0) = p^WLip^m^ = 0)pV., ~ p^ (1.3) 

and its nonrenormalization leads to the fact that the one- loop result w^]^ ~ doesn't 
get the perturbative corrections from gluon exchanges in the higher orders 0|. Nonrenor- 
malization of implies the same for wt'- 

Wl{p\ = 0) = 2wt{p'^, = 0), = 0) = 0. (1.4) 

This nonrenormalization, in difference with longitudinal part, holds on only perturbatively 
and, seemingly, is not directly related to the phenomenon of anomaly. More general non- 
renormalization theorems in perturbative QCD were proved in jl5| . 

It is very interesting to look at the phenomenon of anomaly and perturbative nonrenor- 
malizability of wl,t using the dispersive approach (16i]. In framework of this approach 
anomaly becomes quite simple and represents itself just as an obvious subtraction constant. 
As a result, the two- photon matrix element of the axial- vector current acquires a pole in the 
chiral limit so the anomaly appears as a pure infrared effect. Detailed investigation of the 
one-loop VVA triangle graph within such approach has been performed earlier [l^ . 

The language of dispersion relations allows us to extract some new specific properties of 
the higher order corrections to fermion triangles. Relation ()1.2|) in the context of dispersive 
approach emerges due to the universality of anomaly when it is appearing in the dispersion 
relation in the axial and vector channels. This resembles the mentioned symmetry with 
respect to the permutation u a, observed by A.I. Vainshtein. At the same time, here this 
relation is entirely related to the anomaly phenomenon. 

By direct analytical calculations of relevant formfactors in two-loop approximation taking 
into account symmetry properties of the amplitude we obtained that the full two-loop AVV 
amplitude is equal to zero for arbitrary fermion masses. Non-anomalous Ward identity for 
imaginary parts of formfactors has been proven that provides the correctness of dispersive 
approach usage at least in two-loop approximation. We also make a suggestion that the 
dispersive approach is applicable in any order of perturbation theory. Together with Adler- 
Bardeen theorem immediate consequence of such suggestion is that the Vainshtein's theorem 
is correct for nonzero fermion mass. 

The paper is organized as follows. In Section 2 the dispersive approach to the anomaly 
is briefly described for particular configurations of the external momenta in the standard 
one-loop triangle graph following to [l7'|. The Born approximation to Vainshtein theorem 
in the framework of dispersive approach is interpreted as equality of two expressions for 
axial anomaly and valid in the case of finite fermion mass. In section 3 the generalization 
of dispersion approach to the axial anomaly for higher orders of perturbation is suggested. 
The two-loop radiative corrections to anomalous triangle with arbitrary fermion masses for 
the same kinematical configurations are considered in the context of such an approach. 
We found that all two-loop formfactors are zero, justifying the postulated generalization of 
dispersive approach and Vainshtein's theorem. Section 4 contains some concluding remarks 
and discussion of the higher orders of perturbation theory and non-perturbative effects. 
Appendix contains some details of the two-loop calculations. 
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II. AXIAL ANOMALY AND NONRENORMALIZATION THEOREM IN THE 
FRAMEWORK OF DISPERSIVE APPROACH 

We use the standard tensor representation of the VVA triangle graph amphtude (Fig. El) 
due originally to Rosenberg 0] 

+ k^ea^p^kPp''F5 + p^,eaupak''p''Fe. (2.1) 




FIG. 2: Anomalous triangle diagram. 

Here Fj = Fj{q^\ k"^ , , m^) , j = 1, ...,6 are the Lorentz invariant formfactors. The Bose 
symmetry of the amplitude Ta^^ik^p) = T^ypip, k) is equivalent to 

Fi(/c,p) = -F2{p,k) 

F3{k,p) = -F,{p,k) (2.2) 
F,{k,p) = -F,{p,k). 

The gauge invariance leads to the vector Ward identities k^^T^pu = 0, p^T^pu = which 
in terms of formfactors gives 

Fi = {kp)F^+p'Fi 

F2 = k^F, + {kp)F^. (2.3) 

In one-loop approximation all formfactors Fj may be expressed in Feynman-parametric 
form. For relevant formfactors we have 

F^{k,p) = -—I,^{k,p), F,{k,p) = ^[l2o{k,p) - ho{k,p)], (2.4) 



where 



/•I rl-x x'^y^ 
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Here Nc = 1 for simplicity. These integrals have the following symmetry properties 

Imnik,p) = InmiP, k) (2.6) 

It is useful to observe that ()2.2|) together with ()2.3|) and ()2.6|) implies 

Fe{k,p) = -Fs{k,p) (2.7) 



It is well-known from the classic Adler's paper p that the accurate calculation of loop- 
momentum integrals at the one-loop level leads to the anomalous axial- vector Ward identity 

q"To^f.u{k,p) = 2mT^,{k,p) + (2.8) 

The pseudotensor T^^ may be written as 

T^uik,p) = Gik,p)e^,p^k''p''. (2.9) 
In terms of formfactors ()2.8j) reads 

F2-F1 = 2mG + -^. (2.10) 

According to the Adler-Bardeen theorem the axial anomaly occurs only at one-loop level. 
So, relation (|2.1(Jj) for full formfactors remains the same 

pitot) _ pilot) ^ ^^Q[tot) ^ 1 _ (2.11) 

In the framework of dispersive approach we deal with the imaginary parts of relevant 
formfactors. We start with the kinematical configuration of external momenta fc^ = 0, 7^ 
for particular case where > 4m^ andp^ < 4m^. Then there is a cut for G (4m^, 00) while 
there is no such singularity with respect to variable p^ . For the formfactors Fj, j = 1, ...,6 
and G in any order of perturbation theory one may write unsubtracted dispersion relations 
with respect to 

.2. ^2 _2^ _ \ r ^4/(t^pV«2) 



FM'-,p\rn') = - / ^y^; ' dt, (2.12) 



^,222, i r BHt;p',rn') ^^ 
G{q ■,p ,m) = - I dt 



^2. „2 „„2^ _ 



We use the notations Fj{q'^; p"^ ,m^) = Fj{q^] k"^ = 0,p^,m^) and G{q'^]p^,m^) = G{q^] k"^ = 
0,p^,m^) below in the current section. The and are the corresponding imaginary 
parts, implying the cut with respect to variable = t, for example 

4A. 2.^2 2x _ F,(g^ + ^£;p^m^) - F,{q^ - ie-p\w?) 
/ij [q , p , m ) — — 

The imaginary parts of the relevant formfactors satisfy non-anomalous Ward identities 
because they do not contain the linear divergences in the momentum integrals 

{p^ -t)A^{t]p^,m^)-p^A^{t]p^,m^) = 2m5^(^;p^ m^) (2.13) 
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Using (ITT^ and (ITT!^ one gets finally 

1 f°° 

F2{q'^;p^,m'^) - Fi{q^;p'^,m'^) -2mG{q^;p'^,m^) = - A^{t-p^,m^)dt (2.14) 

Comparing with ()2.1()j) and taking into account ()2.7|1 we find that the occurrence of the axial 
anomaly at one-loop level is equivalent to a "sum rule" [l7l | 



OO -1 

A^{t-p\m^)dt = - (2.15) 

It is easy to evaluate the by taking the imaginary part of the corresponding integral 
in (I23D Q 

l + R\ „ 4m2^^/^ 



Ai{q', p\ m-) = l^j^—Tyl [~P'R + 2m^ In ^ J , = ^1 - ^ j (2.16) 

By integration of this expression over q^ one can check the relation ()2.15|) . 

In the preceding discussion we have employed dispersion relations with respect to variable 
q^ as these were appropriate for considered kinematical region. Let us now discuss another 
version of dispersive calculation of the anomaly using the cuts with respect to p^. For this 
purpose we will consider another kinematical region p^ > Am?, q^ < Am?. Writing now 
instead of (j2.12j) unsubtracted dispersion relations with respect to variable p^ = t for Fj and 
G, we obtain 

F2{q^]p^,m'^) - Fi{q^]p\m^) - 2mG{q^]p\m^) = 
= - / [Al{q';t,m')-Aliq';t,m')]dt, 

The AJ and are the corresponding imaginary parts, implying the cut with respect to 
variable p^ = t 

\V I 2. 2 2\ 



A]{q';p\m') 



Fj{q^]p^ + i£,iT?) — Fj{q^]p'^ — ie.m^) 
2i 



Thus, to recover the standard one-loop anomaly ()2.1U|) taking into account ()2.7|) one has to 
show that 

/ [Al{q^-t,m^)-Al{q^-t,m')]dt = - (2.17) 

JAm^ 27r 

for an arbitrary m and for any considered value of q^. 

A straightforward calculation at one-loop level using ()2.4j) gives a result [l 

\V|2.^2 2^ _ 1 1 /2c 0»^2 ^ 1 + 5* 



Al{q'-y.m^) = ^J^^[PS-^^^^T^)^ (2-18) 

A:{q';p',m') = ^^^5, S= (l-^y (2.19) 
2TT p^ — q^ \ p^ J 
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It was observed [17[ that the integrands A^{q^]p^ ^m?) — A^(g^; j9^, m^) occurring in the 
sum rule ()2.17p at one loop level are equal to the expression for A^(g^;p^,m^) from sum 
rule ()2.15|) with and being interchanged. As a result we have 

AX{q'-p\m') = A^{p'-q\n{') + Al{q^-p\m^). (2.20) 

Let us write the unsubtracted dispersion relations with respect to = t of the both sides 
of (IT^ : 



1 p Al{q^-t,m')^^_ 1 A^{t-q\m^)^^ , 1 Al{q^-t.m 



Jim^ t-P^ TT 74„2 t - p^ VT ^4^2 t - p^ 

One can immediately get from this expression that 

F^{q^-p\ m^) = Fs{p^; q^; m^) + Fs{q^;p^; m^). (2.21) 

In the case with one soft k ^ Q and one virtual photons 7^ we have q^ = (k+p)^ ~ p'^. 
For this kinematics one can obtain the expressions for the longitudinal and transversal parts 
of amplitude Tau in terms of Rosenberg's formfactors 

Wiip^^Tn^) = Fi{p^;p^,m^), 

WTip^m^) = F^{p'-p\w?)-F^{p'-p\w?). (2.22) 

The relation ()1.2|) between wt and in terms of formfactors 

F^{p'-p\m') = 2F,{p^;p^;m^) (2.23) 

immediately follows from ()2.2H1 with = p'^. 

In the framework of Vainshtein's approach the axial anomaly is expressed only through 
the longitudinal part of triangle wl in the chiral limit (jl.3p . Within the dispersive approach 
we have two dispersive relations for axial anomaly ()2.15p . ()2.17p including imaginary parts 
of both structures Wl and Wt for arbitrary mass. 



III. CALCULATION OF TWO-LOOP AXIAL ANOMALY AND CHECK OF DIS- 
PERSIVE APPROACH AND VAINSHTEIN'S THEOREM 

We propose the generalization of dispersion approach to the axial anomaly for any order 
of perturbation theory. We suggest that the imaginary parts of the relevant formfactors 
satisfy non-anomalous Ward identities for higher order of perturbation theory as well. This 
implies that anomaly will be also given by corresponding finite subtractions. 

To check that we will calculate the triangle diagram in two-loop approximation. The re- 
sults for QED and QCD corrections differ only by the obvious colour factor. We consider the 
full amplitude of anomalous triangles Tailu°°^^ with all possible types of radiative corrections 
shown at Fig. El 

At first we construct four possible scalars from amplitude T^^y (j2.ip in the particular 
kinematics A;^ = 0, 7^ 0: 

51 = T^^^k'^e^^'^'^k^p, = -2{kpfF^, 

52 = T^^uP'^e^^'^'k^p, = 2{kpfFi, 

53 = T,^.£"^"''A;, = 6{kp)F2 + 2{kpfF,, 

^4 = T,/..£°^>„ = Q{kp)Fi + 2(^2^5 - 2{kpfF3 + Qp'F^. (3.1) 
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FIG. 3: Two-loop triangle diagrams. 



Together with vector Ward identities for formfactors in the same kinematics 

Fi = {kp)Fs+p''F^ 

F2 = {kp)Fe. (3.2) 

we have the closed system of equations for all formfactors Fj = Fj{q'^;k^ — 0,p^,m^), j = 
1, 6. Its solution is 

p TP {kp)S2-p'{kp)Ss-3p'S, 

^1 — NO ! ^2 — — Trrr-TTTj ^3 — 
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2{kpy' 2{kpf' ^ 2{kp) 

_ {kp)S, + 35i _ S,{kp)-2S2-p^S:, _ 
' 2{kpf ' ' 2(A;p)3 ' ^ 2(fcp)3 ^ ^ 

In the framework of dispersive approach we are interested in calculation of imaginary 
parts of corresponding formfactors that is , p^,m^), so let's calculate the 

imaginary parts of scalars Sn = Sn{q^', k"^ = 0,P^, m'^), n = 1, 4, for example, with respect 
to q'^ 

A^(g2; k^ = 0,p2, m^) = hSniq'' + iS] k^ = 0,^^, m^) _ 5„(g2 _ ^2 ^ q^^2^ ^2)|_ 

We use the Pauli-Villars regularization with the parameter M. After integration over 
loop momenta we have obtained the expressions for rn?) in the form 

of Feynman-parametric integrals (see Appendix). We drop all the terms with the imaginary 
parts with respect to q^ which do not survive in the limit M — > 00. 

The next step is to calculate the imaginary part of Feynman-parametric integrals and 
to get the expressions for An*"^ A;^ = 0,p^,m^) in explicit form before taking off 

the regularization. We get the imaginary part of each integrand and then we integrate the 
(5-f unctions with complicated arguments analytically step by step. After second or third 
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integration all Feynman-parametric integrals in A„ {q'^; k'^ = 0,p^, m^) turned to zero, 

so 

A^(^-'°"P\q^;k^ = 0,p',m'')=0, ri = l,...,4. (3.4) 

According to ()3.3p we see that the imaginary parts of all formfactors with respect to q'^ 
turned to zero in the considering kinematics 

Af-^-'\q';e = 0,p^m')^0, j = l,...,6. (3.5) 

As a result the anomaly sum rule ()2.15|1 preserves one-loop form also at two-loop level. 

The investigation of the Vainshtein's theorem requires the discontinuities of formfactors 
with the respect to variable (or k"^). Unsubtracted dispersion relations (j2.12|) guarantee 

that all two-loop formfactors are equal to zero 

F^'-^'^'^P\q^;k' = 0,p^,m')^0, j = l,...,6. 

So, the full two-loop AVV amplitude in the considering kinematics is equal to zero. Conse- 
quently, the discontinuities with respect to (or k"^) are also zero 

A^^'-^""'\q'-e = 0,p\m')^0, j = l,...,6, 

and the anomaly sum rule ()2.17j) also preserves one-loop form. This immediately leads to 
validity of Vainshtein's theorem with finite mass at two loops 

wl{p^ -iTn^) = '2,wt{p^ -iTTi?')-, w^^j!'°°^\p'^ ,771^) = 0. (3.6) 

For completeness we also explicitly checked the correctness of non-anomalous Ward iden- 
tities for imaginary parts. To do so we have considered the imaginary part of pseudoscalar 
formfactor 5^(2-«oop)(^^2. ^2 

, p^,m^) in the particular kinematics. By direct diagrammatic 
calculations of the two-loop amplitude j'^^"'""*') -^^e found that 

5^(2-'°°P)(g2;A;2 = 0,p2,m2)=0. 

So, taking into account ()3.5|1 . the relation ()2.13|1 does not obtain the perturbative corrections 
in two-loop approximation. 



IV. DISCUSSION 

In our work the axial anomaly and Vainshtein's nonrenormalization theorem are consid- 
ered in the framework of dispersive approach. We found that all two-loop contributions to 
formfactors Fj{q'^;k'^ = 0,j9^,m^) and G{q^\k'^ = 0,p^,m^) are equal to zero for arbitrary 
fermion mass. It allows us to prove the suggested generalization of dispersive approach 
to axial anomaly and to expand the Vainshtein's nonrenormalization theorem for arbitrary 
fermion masses in the triangle loop. 

Although we thus proved these properties only at two-loop level, they are likely to be 
valid at any order of perturbation theory. Indeed, the validity of non-anomalous Ward 
identities for imaginary parts is due to the absence of the linear divergencies and should 
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hold at all loops. This, together with Adler-Bardeen theorem, would result in validity of 
anomaly sum rules. As soon as the corrections to their integrands are zero in chiral limit 
due to Vainshtein's theorem, it is rather hard to imagine the function which is non-zero for 
finite mass case, while its integral is still zero. In turn, all other functions should be also 
zero due to gauge invariance. 

At the same time, the further studies of dispersive approach at higher orders and (es- 
pecially) beyond perturbative theory are of most interest. One should note here the recent 
calculations in the framework of instanton model |18| leading to the exponential, rather than 
power corrections to Vainshtein's theorem. 

It seems that it is the nonlocality of this model, rather than instanton specifics, that pro- 
vides this exponential behavior. In fact, it is analogous to the exponential falloff of transverse 
momentum dependent parton distributions whose coordinate description ^] bears similar- 
ity to the vacuum non-local condensates. One should also recall another observation and 
suggestion of namely, that the local vacuum condensates do not match the dispersive 
description of axial anomaly and that the non-local condensates may improve the situation. 

We are grateful to A.E. Dorokhov, J. Horejsi, M. Knecht, J. Stern, O.V. Tarasov and 
A.I. Vainshtein for useful discussions and comments. We are also indebted to anonymous 
Referee for the valuable comments which allowed to eliminate the inconsistencies in the first 
version of the proof which, however, did not affect the final result. This work is partially 
supported by grant RFBR 03-02-16816. 

Note added 

After finishing of our work we learned about the preprint of F. Jegerlehner and O. V. 
Tarasov j2C?] whose conclusions are similar to ours in the chiral limit for arbitrary external 
momenta. 

V. APPENDIX 

The imaginary parts of invariant scalars (j3.1|) are given by the following Feynman- 
parametric integrals 
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1 1-a 1 1-x 



Af^^~^"°^^ Disc g2 jda j dbjdx j dy [4:{2{kp)Xi - p'^X2){x + 2y - 1)] + 





1 1-a 1 1-x 



+ j da j db J dx j dy [{2>2y{kpf + l&xa{kp)p^ + l&h{kpf - 2>2m'^y{kp) + 



+ 8{kp)p'^ — Sm^yp^ — \Q{kpY + \Qx{kpY + 32ya{kp)p'^ — 8x{kp)p'^ — 32yb{kp)'^ — 

— 16xb{kp)^ — 16m'^x{kp) — \Qy{k:p)p^ — lQa{kp)p^)Xi + (4yp'^ — 8{kp)p^ — IGyap'^ — 

— Sxap"^ — 8m'^x{kp) + Sap"^ — lQm?p^ + 2Am^yp^ + IQni^xp^ + 8y{kp)p^)X2 + {Ap^ — 

— Sx{kp)p^ — Axp^ — Sm^p^ — Ayp"^ + Sm^yp^ + Sm^xp"^ + S{kp)p^ — Sy{kp)p^)X^ + 

1 1— a 1—a—b 1 1—x 

+ J da J db J dcj dx J dy [-16(A;p)^m^(2m^y - 2m^ - p^y + p^ + 2{kp))Xl + 







+ ?>2m'y{kpyX^\ + 

1 1-a 1 1-x 

+ j da j db j dx J dy [Arn^{-^y{kpf - Axa{kp)p^ - Ab{kpf + Aw?y{kp) - 2{kp)p^ + 



+ 2m^yp^ + A{kpY - Ax{kpf - 8ya{kp)p^ + 2x{kp)p^ + 8yb{kpf + Axb{kpf + 

+ 2m^x{kp) + Ay{kp)p^ + Aa{kp)p^ + 2m^{kp))Xi + Am^{-yp'^ + 2{kp)p^ + %aj9^ + 

+ 2xap'^ + 2m'^x{kp) — 2ap'^ + Sm'^p'^ — Am^yp^ — Sm'^xp'^ — 2y{kp)p'^)X2 + 4m^(— p*^ + 

+ 2x{kp)p^ + xp'^ + 2m^p^ + yp'^ — 2m^yp^ — 2m^xp^ — 2{kp)p^ + 2y(A;p)p^)X3] 
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1 l-a 1 l-x 



^A(2-loop) ^ ^ .^^ 



J da J db J dx J dy [-4{2{kp)Xi - p'^X2){x + 2y - 1)]+ 





1 l-a 1 l-x 



+ j da j db I dx j dy [{-32y{kpf + 16{kpf - lQxa{kp)p^ - S{kp)p^ + 
b 

+ lQm^y{kp) + lQ)m^x{kp) + 8x{kp)p'^ — 16x{kpy + lQxb{kpy — 166(/cp)^ + 16y{kp)p'^ + 
+ ?>2yb{kpf - ?>2ya{kp)p^ - Sm^p'^ - 16m^(A;p) + 16a{kp)p^)Xi + {-24y{kp)p^ + 
+ 16yap^ + 9>{kp)p^ - Sm'^xp'^ - Sap^ + 4xp^ - Sm^yp^ - 12yp^ + Sxap^)X2 + (- V + 
+ 4xp^ + 9>y{kp)p^ - Sm'^yp^ + Ayp'^ - Sm^xp^ + 8mV + 9>x{kp)p^ - 8(kp)p^)Xs - 
- Sm^x{kp)Xe] + 

1 1— a 1— a— 6 1 l-x 

J da J db J dc j dx J dy[-l&{kpfm^{-2{kp)y-2m? + 2{kp)^2m?y^p^^ 



+ 







+ xp^)Xl - 16{kp)W{2m^ - 4m^y)Xi] + 

1 l-a 1 l-x 

+ J da J db J dx J dy [-Am^{A{kpf - 4:x{kpf - 8y{kpf + Axb^kpf - 2{kp)p^ + 

DODO 

+ 2x{kp)p^ + 8yb{kpY - Axa{kp)p^ - Sya{kp)p^ + Ay{kp)p^ + 4:a{kp)p^ + 2m^x{kp) - 

- 2mV - 2m^{kp) - 4fe(A;p)^)Xi - 4m^(-2ap^ + 4yap^ + 2{kp)p'^ - 6y{kp)p'^ - 

— m^xp^ — w^p^ + xp^ — 3yp^ + 2xap^)X2 — A.rri^[—p^ + 2x{kp)p^ + xp^ + 2m^p^ + 

+ yp^ ~ 2m^yp^ — 2m?xp^ — 2{kp)p^ + 2y{kp)p^)X^ + 8rn^x{kp)XQ\ 
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1 1-a 1 1-x 



^A(2-loop) ^ ^^^^ 



J da J db J dx J dy [-8{Xi + X2){x + 2y - 1)]+ 





I— a 1 1—x 



+ J da j db I dx j dy [{-32{kp)y + 32yb{kp) - IQxap^ - 32yap^ + 16p^y - 



- 16x{kp) - + 16m^x + 8xp^ + 16{kp) + 32m^y + 16m^ - 16b{kp) + 16ap^ + 
+ 16xb{kp))Xi + {-32{kp) + Sxp"^ + AOp^y + Idrn^x + 16m^y - 16p^ + 16x{kp) + 
+ 16ap^ + AS{kp)y - 32yap^ - 16xap^)X2 + {Sp'^y + 16m^x + 16m^y - 16{kp)y + 
+ 8p^ -16m^ + 16{kp))X3] + 

1 1— a 1—a—b 1 l—x 

J da J db J dc j dx j dy [16m'^{kp){-2{kp)y - im^ + im^y + 4:{kp) + 2p^ - 



+ 





2p'^y + xp'^)Xl - 64m'^y{kp)Xl] + 

1 l-a 1 1-x 



+ J da J db j dx j dy [—8m'^{—2x{kp) — 4:{kp)y + Ayb{kp) — Ayap^ + 2p^y — p^ -\- 
b b 

+ m^x + xp^ + 2{kp) + 2w?y + 3m^ - 2b{kp) + 2a/ + 2xb{kp) - 2xap^)Xi - 
- Sm'^{5p'^y - 2p^ + rn^x + xp^ + - Ayap^ - A{kp) + 2ap^ + 6{kp)y - 2xap^ + 

+ 2x{kp))X2 - Sm^{2m^x - 2m^ + 2m^y-p^y + 2{kp) - 2(A;p)y + p^jXg] 
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1 1-a 1 1-x 



^Ai2-loop) ^ ^ .^^ 



J da J db J dx J dy [8{-l + 2y + x)Xi + 8{-l + y + x)Xe]+ 





1 1-a 1 1-x 



+ j da j db I dx j dy [(24/ + lQ{kp) - lQx{kp) - 16m^ + m{kp) - 24^ - 
b 

- 32m^y - 32yb{kp) + 32yap^ - 16ap^ + Idxap'^ - ?>2p^y - 16m^x - 16xb{kp))Xi - 

- IGyp^Xa + {Sp^y - l&m^x - 16(A;p) + 16m^ + l&{kp)y - - lQm'^y)X^ + 

+ {-l%{kp) - l&x{kp) + lQyb{kp) + l&xb{kp) - 16p^y + 16p^ + ?>2{kp) - 32{kp)y - 

- 16m'^x)XQ] + 

. j ^ J\ / (8(-e™VM/ . 8.^M/ . - 



— 8m'^y{kp) — 2x{kp)p'^ — 8rn?{kpY + ^qj^kpYp^ — p^ + 2m^xp^ + 8m^y{kpY + 
+ 8m^{kp) + 8m^yp'^ — A{kp)p'^ — A^kpYp^ + yp^ — xp^ + Qrn^p'^ — Am^x{kp)p^ — 

— — Qm^yp'^)Xl + 8{—Am^y{kp)'^ + 2m^{kp)p^ + 2m'^yp'^ — 8m'^{kp) — 

— 8m^yp^ — Am^x{kp)p^ — Am^x{kp)^ + I6m'^y{kp) + 2m^y{kp)p^ + Am^{kp)'^)X^] + 



1— a 1 1— x 



+ y rfa J db J dx J dy [8m^{-2{kp) + 2x{kp) + 3m^ - 26(/ep) - + 2m^y + 
0000 

+ Ayb{kp) - Ayap^ + 2ap^ - 2xap^ + Ap^y + m^x + 2xb{kp) + 3xp^)Xi + IGm^yp^Xa + 
+ 8m^(2m^x - 2m^ + 2m^y - p^y + 2(A;p) - 2{kp)y+p^)Xz + 8m^(-2x6(A;p) + 26(A;p) + 

+ 2p^y + — m^y + 2x{kp) + m^x — 4(A;p) — 2yb{kp) — 2p^ + 4(/cp)|/)X6] 
where 

Xi = 1/ [—rn^ya — m^yb + m^y — rr? + 2ab{kp) + p^a — p'^a^); 

X2 = 1/ (rn^xa + m^a;6 — m^x + rn^ya + — m^y — m^a — m^b — 2xya{kp) — 

— 2xyb{kp) + 2xy{kp) — p^xa — p^xb + p^x + p^x^a + p^a;^6 — p^x^ + p^a — p'^a^); 
Xs = 1/ {m^xa + m^xfe — m^x + m^ya + m^|/6 — w?y — w?a — m% — 2xya{kp) — 

— 2xyb{kp) + 2xy{kp) + 2ab{kp) — p^xa — p^xb + p^x + p^x^a + p^a;^6 — p^x^ + 

I 2 2 2\ 

+ p a — p a ); 

= 1/ {m^xa + 777,^x6 + m^xc — m^x + m^ya + m^yb + m^yc — m^y — m^a — m% — 

— m^c — 2xya{kp) — 2xyb{kp) — 2xyc{kp) + 2xy{kp) + 2ac{kp) — p^xa — p^xb — p^xc + 

,2 ,22 ,2 2i, ,22 22,2 2 2\ 

-\-px-\-pxa-\-pxb-\-pxc — px + p a — p a ); 
Xq — l/{—m^ya — rn^yb — m^yc + m^y — rr? + 2ac{kp) -\-p^a — p^a^)] 
Xq — 1 / (2xy{kp) — 2xyb{kp) — 2xya{kp) — m% — m^a — m^y + rn^yb + m^ya — rr?x + 
+ vr?xb + rr?xa + p^x — p^xb — p^xa — p^x^ + p^x^ + p^x^a). 
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